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We construct static and spherically symmetric black hole 
solutions in the Einstein-Euler-Heisenberg (EEH) system 
which is considered as an effective action of a superstring the- 
ory. We considered electrically charged, magnetically charged 
and dyon solutions. We can solve analytically for the magnet- 
ically charged case. We find that they have some remarkable 
properties about causality and black hole thermodynamics de- 
pending on the coupling constant of the EH theory a and b, 
though they have central singularity as in the Schwarzschild 
black hole. We restrict a > because it is natural if we think 
EH theory as a low energy limit of the Born-Infeld (BI) the- 
ory. 

(i) For the magnetically charged case, whether or not the 
extreme solution exists depends on the critical parameter 
a = acrit- For a < acrit, there is an extreme solution as in the 
Reissner-Nortstrom (RN) solution. The main difference from 
the RN solution is that there appear solutions below the hori- 
zon radius of the extreme solution and they exist till rn — > 0. 
Moreover, for a > acrit, there is not an extreme solution. For 
arbitrary a, the temperature diverges in the rn ^ limit. 

(ii) For the electrically charged case, the inner horizon ap- 
pears under some critical mass Mo and the extreme solution 
always exists. The lower limit of the horizon radius decreases 
when the coupling constant a increases. 

(iii) For the dyon case, we expect variety of properties be- 
cause of the term bief^^paF^'^ F'"')'^ which is peculiar to the 
EH theory. But their properties are mainly decided by the 
combination of the parameters a + 86. We show that solu- 
tions have similar properties to the magnetically charged case 
in the rn Q limit for a + 86 < 0. For a + 86 > 0, it depends 
on the parameters a, 6. 

04.70.-S, 04.40.-b, 95.30.Tg. 97.60.Lf. 



I. INTRODUCTION 

Recently, much attention has been payed to BI type 
of actions after its recognition as an effective theory of 
superstring theory j^]. Moreover, since they describe the 
action of the Brane, their importance has been increas- 
ing Q. In this context, there are some studies that in- 
vestigated black hole solutions in the Einstein-BI type 
actions [p|. Actually a new non-linear electromagnetism 
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was proposed, which produces a nonsingular exact black 
hole solution satisfying the weak energy condition 
and have distinct properties from Bardeen black holes 
But there remain subjects which should be manifested 
such as thermodynamical properties. Here, we concen- 
trate on the EH action which first was proposed in 1936 
. Though not so much attention has been payed it com- 
pared with the BI action, EH action well approximate 
the supersymmetric system of minimally coupled spin- 
1/2, -0 particles for appropriate parameters From 
the experimental aspects, this is a more accurate classi- 
cal approximation of QED than Maxwell's theory when 
fields have high intensity ^ . 

We investigate the black hole solutions in the EEH sys- 
tem from following aspects ; (i) The electric-magnetic du- 
ality (ii) The black hole thermodynamics (iii) The causal- 
ity and stability of the black holes. As for the electric- 
magnetic duality, it is already pointed out that though BI 
action preserves this, EH action breaks it at the higher 
order of the electromagnetic field |l^ . Thus black hole 
solution with electric charge or magnetic charge will have 
clear differences which should be clarified. Moreover, the 
dyon solution may have specific properties which can 
not be seen in the electrically charged or magnetically 
charged cases. The thermodynamical properties of black 
holes are one of the main topics in superstring theory af- 
ter the discoveries of the microscopic origin of the black 
hole entropy |ll) and the holographic principle ||l^ . It is 
worth noting that BI type action also plays an important 
role in AdS/CFT correspondence |l^. Causality for the 
black hole in the EEH system has already been inves- 
tigated by Oliveira. But this is restricted to the black 
hole with electric charge, and the physical implications 
such as the stability of the black hole are not discussed. 
The stability of the black hole can be interpreted from 
thermodynamical properties established in and we 
refer to this by calling it a turning point method. Using 
thermodynamical variables, we can easily apply catas- 
trophe theory to hairly black holes and this is consistent 
with linear perturbation analysis and the turning point 
method [T5| ]. We are interested in the thermodynamical 
properties of black holes in this system and the relation 
between Causality and the stability of black holes. 

This paper is organized as follows. In Sec. II, we intro- 
duce basic Ansatze and the field equations in the EEH 
system. In Sec. Ill, we investigate the thermodynami- 
cal properties of black holes with electric charge or mag- 
netic charge. In Sec. IV, we investigate those of dyonic 
ones. In Sec. V, we summerize the results and com- 
ment on future work. Throughout this paper we use units 
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c = h = 1. Notations and definitions such as Christof- 
fel symbols and curvature follow Misner-Thorne- Wheeler 



II. BASIC EQUATIONS 



We take the following EEH action 
1 



S = d xy/—g 



IBtt V G 



(1) 



where R is the scalar curvature, P 
e 

unit tensor, which yields 



fJ,V7 Q 

fivpcrF^^ F^'^ and e^^^"^ is a completely antisymmetric 



= -4!. 



(2) 



In 1^, this corresponds to the weak field approxima- 
tion and the coupling constants are written as a = 
he^/imOir'^m'^), b = 7he'^ / {lUOn'^m'^), where h, e and 
m are the Planck constant, electron charge and electron 
mass, respectively. From the present point of view, they 
should be related to the inverse string tension a' which re- 
strict a > because of its correspondence to the BI action 
in the low energy limit. As has been pointed out in [Q, 
to construct gravitational counterpart of the BI action 
involves much difficulty. Here, we adopt the Einstein- 
Hilbert action as a first approximation. We can derive 
Einstein equations as 
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G^. = ^9^A-P + aP' + bQ') + 2F^xF^ 
-4aP{F^xF^) - 8bQ{e^cv^F^''P^)- 



(3) 



We consider the metric of static and spherically symmet- 
ric, 



(4) 



where f{r) = 1 — 2Gm{r)/r. We introduce the gauge 
potential A^, as 



A„ 



(A(r),O,O,Q„cos0). 



Then, the Einstein equations are 
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Fdy, 



(5) 



(6) 
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S'f = -Fe - F„ - Fdy, (7) 



r 

where ' represents d/dr. We used the abbreviations as 

2Sr Ai\1 I e„„4(5/ \t\\ 



T? _ Qm r)„ Qm 



Fiy={mb-^a)^\^f^. 



(8) 
(9) 

(10) 



Subtracting Eq. (||) from Eq. yields d5 / dr = 0. We 
require asymptotically flatness for the solution. 



A{r) 



-, 5{r) —> 0, Gra{r) — s- constant, (11) 



as r — > oo. Thus we obtain 5 — Q. So we have only one 
independent Einstein equation as 



Gm' = —{F,+ F„, + Fdy) 

The field equation is 

4ar^A'f + A'z{r) = Q„ 

where 

z{r)=r^ -4{a + 8b)^. 



(12) 



(13) 



(14) 



This is third order algebraic equation for A' except for 
Qe — 0- For regularity at the horizon th, we require 

Gm{rH) = ^TH, A{rH) < oo. (15) 



III. BLACK HOLE SOLUTIONS WITH 
ELECTRIC OR MAGNETIC CHARGE 

In this section, we show the properties of black hole so- 
lutions with magnetic or electric charge. First, we point 
out that the zeroth and the first law of black hole ther- 
modynamics can be applicable even for non-linear matter 
terms which violate dominant energy condition though 
Smarr's formula can not p7[ |. 



A. Magnetically charged case 

In the case Qe = 0, we can solve equations analytic ally . 
In this case, there remains only Fm part in Eq. (|l^). 
Note that Gm' can be negative which makes an intrinsic 
difference from the RN solution. We can integrate Eq. 



Gm = GM - +a-^, 
2r 5r° 



(16) 



where M is the gravitational mass of the black hole. 
Thus, the horizon radius rn must satisfy 

hirn) ^r%~ 2GMr% + Qlrj, - ^aQt = 0. (17) 

5 

Since h{0) < and h(oo) oo, the solution which satis- 
fies h{rH)_^ for rn > always exists. 
From (^), 
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5GMrH + 2Qi 



(18) 



So we can classify the number of the horizon as foUows. 
For L = {bGMf - 24Q^ > 0, if h{[bGM + VI]/6) < 
and h{[bGM — \/L]/6) > 0, there are three positive 
solutions, which means that there are one outer horizon 
and two inner horizons. If /i([5GM + •\/L]/6) = or 
h{\bGM - yi]/6) = 0, there are two horizons. In other 
cases, there is only one horizon. 
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FIG. 1. M-m relation for Q,n/lp = 1 and a/Ql^ = 0, 0.01, 
2/27, 0.1, 1. The points A correspond to the extreme solu- 
tions. The lines between AtoB,BtoC correspond to outer 
inner horizon and inner inner horizon, respectively. We can 
see that acrit = 2Q^/27 divides the properties qualitatively. 
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FIG. 2. M-l/T relation for the same parameters in Fig. 
1. For a < acrit, there is an extreme solution where the 
temperature is zero. For arbitrary a, the temperature diverges 
in the rj/ ^ limit. 



We also evaluate that what condition would be re- 
quired to exist an extreme solution. Gm' = 1/2 leads 
to 



Thus 



K{r) = - Qlr^ + 2aQt, 



(19) 



(20) 



So r = y/2/3Qm is a local minimum of K for r > 0. 
K{^/2/3Qm) = leads to a = {2/27)Q^^ = acrit, which 
means that there is not an extreme solution for a > acrit ■ 
We first show the relation between the gravitational 
mass M and the horizon rn for Qm/lp = 1 and ajQl^ = 
0, 0.01, 2/27, 0.1, 1 (Fig. 1). is the Planck length. For 
'^IQ'rn — 0, 0.01, 2/27, there is an extreme solution (the 
point A). The lines between A to B, B to C correspond 
to the outer inner horizon and the inner inner horizon, re- 
spectively. Note that below the point C, there are black 
hole solutions again. For a/Q^ = 0.1, 1, there is not an 
extreme solution as we showed. In all cases, M ~oo 
for vh 0. We also show the gravitational mass M and 
the inverse temperature 1/T relation in Fig. 2. Above 
the mass corresponding to the point A, it is similar to 
the RN's qualitatively. But below the mass correspond- 
ing to the point C is quite different from RN's. The 
temperature is finite but nonzero at the point G. The 
curve from G to B means that if we regard an inner hori- 
zon as an event horizon, the 'temperature' goes to zero 
when rn approaches the point B. If we apply the turn- 
ing point method in this case, the line below the point 
C would be unstable. So we can regard this as unphys- 
ical. But for a/Q^^ = 0.1, 1, this method suggest that 
there is no stability change if we think in the isolated 
system though thermodynamical properties are different 
in these two cases. The specific heat of the black hole 
never changes for a/Q^^ — 1, while it changes twice at 
the points D and E for a/Q^ = 0.1. In all cases, the 
temperature diverges for rn ~* 0. It is reasonable that 
higher order curvature terms would change the results 
in this region. But even if we believe that this system 
describes black hole solutions correctly, it is difficult to 
observe the negative mass black holes since it will evap- 
orate very quickly. 



B. Electrically charged case 



In this case, from Eq. (|1J 

iar^A'f +A'r^ = Qc 



(21) 



which has only one real solution and two imaginary so- 
lutions. The real solution is 



A\r) = 



-2 •3i/3r2 + 21/3^2/3 



62/3^^^1/3 

We used an abbreviation as 

B = gV^Qe + y/l2r^ + SlxQl, 
x{r) = Aar^ . 



(22) 



(23) 
(24) 
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There remains only Fg part in Eq. which shows 

Gm! > 0. This is one of the main difFcrcnce from the 
magneticahy charged case. We show that an extreme 
solution always exists. If we take G'm'(rjj) = 1/2, A'{rH) 
is evaluated from Eq. dia) as 



A'irH) = 



(y- 1)1/2 



2(3a)i/2 ■ 

We introduced a dimensionless valiable y as 




(25) 



(26) 



Substituting Eq. (|2|) into Eq. (p2|) derives 

g{y) ^ (y - if/' + 3{y ~ 1)'/' ■ 
Thus, 



4(3a)i/2 

QeV 



2 (y- 1)1/2 2(3a)i/2 



. (27) 



(28) 



Because of y > 1, the solution of g'{y) — is y = ya = 
1 + 27a/Ql{> 1). We can see g'{y) > for 1 < y < ?/o 
and g'{y) < for yo < y. So if we notice that g{l) — 
and (7 — > — oo (y — > oo), we find there is only one positive 
solution. So we can conclude that there is one extremal 
black hole for a > 0. 

Electrically charged case has already been investigated 
previously |lq ]. The solution can be expressed using the 
hypergeometric function. But we need numerical calcu- 
lation to investigate their detailed properties, particu- 
larly their thermodynamical properties. The inner hori- 
zon only appears for black hole solution M < Mq as he 
showed. Afo is 



3/2 



r(i/4) Q 

2r(3/2) (2a)i/4' 



(29) 



We first show the field distributions of the solutions ((a)r- 
m (b)r-A') for rn/lp = 1, Qe/lp = 1 and a/Ql = 0, 
0.1, 1, 10 in Fig. 3. Because its difference from the 
Maxwell field is particularly large near the event horizon, 
this makes a nontrivial change for a small black hole. We 
can evaluate that for a^^'Qe S> r^, 



A' 



Q 



1/3 



(2^)2/3^1/3 ■ 



(30) 




(b) 




FIG. 3. Field distributions of black holes with electric 
charge for rh/lp ~ 1, Qe/lp = 1 and a/Q^ — 0, 0.1, 1, 10 
((a) r-m (b) r-A'). Because of the difference from Maxwell 
field at small scale, the resulting solution deviates from the 
RN black hole near the horizon. A' monotonically decreases 
as r — » CX3 as is easily shown. 



Differentiating Eq. (g2|) shows A" < 0, so A' mono- 
tonically decreases as r ^ oo. We investigate the M-th 
and M-l/T relations for electrically charged black holes 
in Fig. 4 and 5, respectively. We take a/Ql = 0, 0.1, 1, 
10. The point A corresponds to the extreme solution and 
the curve ^ to i? shows an inner horizon. For finite a, 
there exists an extreme solution as we noted above and 
this approaches to r// — > for a — > oo. Thermodynami- 
cal properties are similar to the case for the RN solution. 
The point D corresponds to the point where the specific 
heat changes and this is not equivalent to the point B. 
So there is no relation between the point which is rele- 
vant to the causality change and the point at which the 
specific heat changes. 
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FIG. 4. M-rH relation for a/Ql = 0, 0.1, 1, 10. The causal- 



ity changes at the point B below which the inner horizon ap- 
pears. Though the lower limit of the horizon decreases as we 
take a large, the extreme solution always exists. 



A'{r) 



(1 ± i^/i)x 
22/331/3^51/3 



(lT^\/3)gi/3 
24/332/3 ■ 



In this case, 



(32) 



(33) 



Note that for 12z' 
is (P). For a + 8b 



' + SlxQl > 0, the only real solution 
> 0, because z can be negative for 



a small r, 12z'^ + SlxQj. < is possible only near the 
horizon. But even in that case, there is only one positive 
solution (^). We should take a positive solution because 
z eventually becomes positive for large values of r. So we 
take (31) in any case. 

We can classify solutions in the ^ limit three 
types as follows. 

(I) If a + 8& = 0, A'{r) approach (|^ for a^^^Qe > r^, 
which is the same as in the electrically charged case. On 
the contrary, Eq. (|l^) approach 



um ~ —a — 7^, 



(34) 




FIG. 5. M-l/T relation for the same parameters in Fig. 
4. The extreme solution always exists where the temperature 
becomes 0. So the M-l/T relation is similar to the one for 
RN black hole. Note that the point where the sign of the 
specific heat changes does not necessarily correspond to the 
point B which suggests that the causality change will not be 
irrelevant to the stability change. 



IV. DYON BLACK HOLE 

As we showed above, the properties of black holes have 
very different aspects, depending on whether it has elec- 
tric charge or magnetic charge. In this section, one of the 
main purposes is to survey how thermodynamical prop- 
erties change when we change the Qm/Qe ratio or the 
coupling constants a, b. In this case, from Eq. ([l3|), the 
three solutions are expressed as 

_2- 31/^2 -)- 2^/3 R2/3 



which has same form as in the magnetically charged case. 
So the characteristic feature of small rn is like that of the 
magnetically charged case. 

For a + 86 7^ 0, we can see its nature if we rewrite ( |3l| ) 

as 



A' = 



Al-Al 



62/37iBl/3(A2 + A1A2 + Al) ' 



where 



Al = 21/3^2/3^ = 2-31/3 



(II) For a + 86 < 0, we can evaluate 



(35) 



(36) 



(37) 



which shows A' cx in the r ^ limit. 

So we can conclude that if (a + 86) < 0, Eq. (|l|) has 
same asymptotically form ( p^ ) in the rn limit as in 
the magnetically charged case. 

(Ill) For a -I- 86 > 0, there exists r = tq below which 
12z'^ + SlQlx < is satisfied. For a while we consider 
r < To case. Then we can evaluate 



Al 



A^ = -24z^ 



(38) 



62/37^51/3 



in the r limit. So we can not 
conclude whether or not solutions in the rn limit exists. 
It depends on a, 6 as we see below. 

Next, we show the field distributions in Fig. 6 (a) r-m 
(b) r-A' for a/Ql = 1, b/Ql = -1, (i.e., a + 8b < 0), 
rn/lp = 1, Qe/lp = 1 and Qm/lp — 10"'', 1. Monotoni- 
cally decrease of A' is broken and m' < region is specific 
for Qm/lp = 1 contrary to the case for Qmjlp = 10^^. 
But they are universal in the ru limit unless Qm ^ 
as is shown above. 
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(a) 




Q /Q=10'" 

m e 



Q /Q=l 

m e 



charged case. But for b/Ql = 0, 0.1, we can see the 
character like the electrically charged case, i.e., solutions 
below an extreme solution do not exist. The curve below 
the points A is a sequence of inner horizons. We also 
investigated those for various Qe/Qm ratio which suggest 
that whether or not solutions in the r/f — > limit exist 
depends only on a, b. Thus if we believe that this system 
is realistic, the coupling constants decide the final fate of 
black holes. 



10 10 10 10 

r/r 

H 




Field distributions of 



dyon black holes for 
-1, Qe/lp = 1 and 



FIG. 

rn/lp = 1, a/Ql = 1, b/Ql = - 
Qm/lp ~ 10~^, 1 ((a) r-m (b) r-A'). As we can see for 
Qm/lp ~ 1 monotonically decreasing of A' is broken and 
m' < region appears which can be seen in the rn ^ 
limit unless Qm 7^ 0. 

We also studied th and 1/T relations in terms of M for 
above three cases. We first show those for a + 86 < in 
Fig. 7 (a) and (b), respectively. We fixed the parameters 
a/Ql = 1, b/Ql = -1, Qe/lp = 1 and Q^/lp = 10-^ 1. 
For Qm/lp = 1, we can easily see specific properties of the 
magnetically charged case though for Qm/lp = 10~* we 
can not. But it is not true. Even for Qm/lp = 10"'*, there 
exist solutions in the r// — > limit where the temperature 
diverges and M —00. They are clear from Fig. 8 which 
is a magnification of Fig. 7 (a) . 

We show corresponding diagrams for a + 86 > in 
Fig. 9 (a) and (b) , respectively. We fixed the parameters 
a/Ql = 1, b/Ql = -0.125 (i.e., a + 86 = 0), -0.1, 0, 0.1, 
Qe/lp = 1 and Qm/lp = 1- For b/Ql = -0.125, -0.1, 
it is almost indistinguishable in this diagram though the 
electric field has different limit for — > in these two 
cases. We can see the character like the magnetically 




FIG. 7. (a) M-rn (b) M-l/T relations for dyon black holes 



for a/Qi = 1, b/Ql 



-1, (i.e., a + 8b < 0) Qe/lp = 1 and 



Qm/lp = 10-*, 1. 

only exist for Q„ 



It seems that solutions in the rjj 
/Ip = 1. But it is not true. 



limit 
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(a) 
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P 

FIG. 8. Magnification of Fig. 7 (a) wfiidi sliows that there 
exist solutions in the rif — > hmit even for Qm/lp = 10^*. 



V. CONCLUSION AND DISCUSSION 

We note our conclusions and future work. We investi- 
gate black hole solutions in the EEH system for electri- 
cally charged, magnetically charged and dyonic solutions. 
They have remarkable thermodynamical properties. 

(i) For the magnetically charged case, the properties of 
the black holes change qualitatively for a — Ocrit- There 
is an extreme solution only for a < acrit- There are 
solutions in the rjj ^ limit for arbitrary a and the 
temperature diverges in this limit. 

(ii) For the electrically charged case which was ana- 
lyzed previously, though the lower limit of the horizon 
becomes small as we take a to be large, the final state of 
the black hole when we consider the evaporating process 
is similar to the RN one. The causality change is also 
already pointed out, i.e., inner horizon appears only for 
M < Mq, as we confirmed it. But the point M = Mq is 
not relevant to the change of the stability if we apply the 
turning point method. 

(iii) As for the dyon case, we showed that there exists 
solutions in the th limit for a-|-86 < and approach 
the magnetially charged case. For a -|- 86 > 0, our results 
suggest that whether or not solutions in the r// — > 
limit exists depends only on a and b not on Qm/Qe ratio 
except for vanishing Qm or Qe- 

We now comment on our future work. Though we 
considered the Einstcin-Hilbert action as a gravitational 
part, it is important to generalize to think higher or- 
der curvature corrections. It may be interesting to think 
about black hole solutions in the action which generalize 
the EH action to preserve supersymmetry . Our solu- 
tions have pathological properties like the negative gravi- 
tational mass. There may exist mechanism which prevent 
such properties as in pO| . An other concern we have is 
to think about black hole solutions including cosmolog- 
ical term, because its importance is recognized both in 
observatioal and in theoretical perspectives. 
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b/Q„ =-0.125 
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FIG. 9. (a) M-VH (b) M-l/T relations for dyon black holes 
for a/Ql = 1, b/Ql = -0.125, -0.1, 0, 0.1, QJlp = 1 and 
Qm/lp = 1. For &/Qe = —0.125, —0.1, these figures are almot 
indistinguishable for these two cases and resemble those for 
the magnetically charged case. But for b/Q^ — 0, 0.1 they 
resemble those for the electrically charged case. 
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